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Abstract 

This paper ex.ends Farrenkopf s analys.s of a stag, e-axis 

data to include the angle output white noise of a rate Integra 1 g & equations. It is shown that only 

Introduction 

There has fopg been an interest in Ka^an fiheHng- *££££ £ 

attitude estimation. 3 4 Farrenkopf found an analytic solution J useful f J preliminary analysis of 

filter combining data from a gyro and an angle senso , assumed a rate gyro model, however, and many 

spacecraft attitude determination systems. Farrenkopf lmphc y js tQ modify Farrenkopf s 

missions employ highly accurate rate-integrating gyros CRI^ resu |ts are significantly more complex 

ISSSSSS ^variance am , simp,e motion ,0 his. 

We begin by discuss, ng the dynamic models fo, ft. ^ 
star hacker measurement. Then we derive the ” Sons am”er,»ed new. and i, is 

propagated to future times. The conclusions are stated at the end of the paper. 

Tlvnamics 


The single-axis spacecraft dynamics are given by 


0 = CD, 


(i) 


where Sis dre reunion angle and mis the hue angular velocity. Because of its own interna, ‘ d “ s 

no, measure 8 exaedy. but instead accumulates it own angle #. The RIG dynanne equahon for dus angle 


0 = cu + f> + n„, 

Gaussian white noise process. The drift rate is assumed to satisfy 


(2) 


where b is the gyro drift rate and n v is a zero mean 

b = n t , 

where is also a aero mean Gauss, an wh„e noise process. These processes ate assumed to obey 

_2 


E{n v (fK(t')} = cr 3 <5(r-r'). 


and 


£K('K(f')} = o, 

where £(.} denotes the expectation value and 5(.) is the Dirac delta function. 


(3) 

(4a) 

(4b) 

(4c) 


Farrenkopf took 0 to be the measured gyro rate and eliminated co between Eqs. (1) and (2). We will not make this 

assumption, since it does not provide a convenient procedure for including angle output white 35 
readout noise or electronic noise) on the RIG output. 6 The procedure we adopt requires a dynamic model for ^ 


LiVJiat*/ oi t cj ic I'*'-' ^ r ■ l u C I or* 

angular rate. Since we are interested in applications that use the gyros for rate measurements rather than Eu . 
dynamic models, we approximate the rate as a pure random walk, 

(b = n w , 

where n w is a third zero mean Gaussian white noise process, uncorrelated with n u and n v , and obeying 

E{nJt)n w (t')} = a 2 w S(t-l'). 


(5) 


(6) 


We will eventually consider the limit of very large a w . 

The dynamic variables are the four components of the state vector 

x = [d,b,<p,(o] r , 


(7) 


which obeys the dynamic equation 


x = Ax +[0 1 n u ,n v ,n w ] r , 


( 8 ) 


with 


A = 


0 0 0 f 
0 0 0 0 
0 10 1 
0 0 0 0 


(9) 


This state vector contains two more components than Farrenkopf s, which only contains 9 and b. With a caret 
denoting expectation values, Eq. (8) gives the state propagation equation 

i(r) = <t>(r-r 0 )x(r 0 ). 


( 10 ) 


where Of / -t 0 ) is the state transition matrix 




TOO t-t 0 

0 10 0 

0 t - 1 0 1 t - 1 0 

0 0 0 1 


(ID 


The state covariance matrix 


/> = £((.r-x)(x-T) T ) 


( 12 ) 


obeys 


P(t) = <D(r - t Q )P(t 0 )<I> T (t - 1 0 ) + f <*>(' " O' Q ^ if ~ f ') - 


(13) 


where 

Qsdiag([0,CT„\crJ,cr;]) 


( 14 ) 


is the process noise spectral density matrix. 



Measurement Models 


Our estimator processes gyro measurements at a time interval r, and star sensor measurements at a multiple of this, 

T = nr. ( 15 ) 

Since we are looking for a steady-state solution, we will only consider a single star tracker measurement followed by 
n gyro measurements. We will require the covariance after processing all these measurements to be identical to the 
covariance at the beginning of the time interval of length T, which is assumed to be immediately before the star 
tracker measurement at time zero. 

The star tracker measurement at time zero is modeled as the spacecraft rotation angle plus a white noise term. 

6 0 =6(0) + n n , (16) 

where n n is the star tracker measurement noise, which is assumed to obey 

E{V>/.} = a l S r (17) 

The sensitivity matrix for the star tracker measurement is the row vector 

H a =[ 1, 0,0,0]. (18) 

At time it, the gyro provides a measurement which is its integrated angle plus a white noise term: 

<p. = 0(;r) + n, f , for i= 1,2, ... ,n, (19) 

where (p i is the measured value and n* is the gyro measurement noise, which is assumed to obey 

E{n‘n]} = c]8 ir (20) 

with <5. denoting the Kronecker delta. The gyro noise and the star tracker noise are assumed to be uncorrelated with 
each other and with the process noise. The sensitivity matrix for the gyro measurement is the row vector 

H g =[0,0, 1,0]. (2D 


These are the models we need to derive our results. 


Gyro Measurement Processing 

Assume that we have performed a star tracker update at time zero, so that x(0+) and P(0+) denote the state 
expectation and covariance immediately after a star tracker update. Then i(/T—) and P(iT-) for i = 1,2, ... , n will 
denote the state expectation and covariance immediately before the gyro update at time it, and x(ir+) and P(i r+) 
will denote these quantities immediately after the gyro update. The steady-state condition is that P(nr+), the 
covariance immediately before the star tracker update at time T, is equal to P( 0-), the covariance immediately before 
the star tracker update at time zero but after the gyro update immediately preceding it. 


Denoting the components of the post-gyro-update covariance by lower case p, 
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Eqs. (11M14) give 
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The standard Kalman filter equations for the gyro measurement update give 1 - 2 

P(ir+) = P(iT-) - P(iT-)H]A~'H t P(iT-), 


( 23 ) 


(24) 


where A, is the scalar 


Ai = H t P(ix-)Hl + a] ^ ^ (25) 

- P„U - 1) + 2x[ Ptf (i - 1) + P^d ~ 1)] + r 2 [ Pbb (i - 1) + 2 Pb Ji - 1) + Pmm O' - D1 + ra 2 v + i x 3 (cr 2 + <r 2 ) + <r 2 . 

Combining Eqs. (22H25) gives the ten scalar equations 


Pge (i) = Pee (i - 1) + 2r Pga> (i - 1) + x 2 Paa) (i - 1) + i^a 2 

- A' 1 {/>«*(« - 1) + X[ Peo (i - 1) + Peb (i - 1) + fV»(* - !)] + V + - !)] + 

Peb (i) = Peb(i - 1) + r Pba>(i ~ U _ A i' ( Pa *(* “ *) + T [Pto(* ” ^ + P®(* ~ *) + P*»(‘ ~ ^ 

+ t 2 [Pm,(‘ " 1) + A*u(» " 1)1 + iT 3 cr 2 }{p w (i - 1) + x [ Pbb (i - 1 ) + PtmO ~ 1)] + |T 2 ^ 2 }. 

Pbb (0 = Pbb 0 ~ 1) + T<7 2 ~ A’ 1 [p M (i ~ 1) + x [ Pbb O' - 1) + Ptc U ~ 1)1 + i T 2 ^ 2 } 2 > 

p H { i) = <T 2 A- 1 [ Pat {i - 1) + x[ Peu (i - 1) + P ®( i - 1) + P+.V - D] + T2 tP»»(i - D ■ + /»-(*' - Dr + i tW ~ ] ’ 

p*0') = A i' (PmO' - 1) + *lp»0 - 1) + Pbc, O' - !)] + I t2£T « 

p*(0 = <^ 2 a; 1 (a, -o- 2 ), 

PeJO = Pm, O' ~ 1) + ^ P m 0 ~ 0 + + “ A ^‘ + " !) + P *“ ( ‘ " 1)1 + X 

{/>„,(/ - 1) + T[ Pew d-l) + Pgb (i - 1) + P*,(« - 1)] + r 2 [p*u(« - !) + P -»(*■ - l » + 

PbotU) = Pbw U - 1) - A ,' {p*,(« - 1) + r [Pbw^ - !) + Po , «(* " !)] + i r2<r »} x 

{p*(i - 1) + r[p«0' - 1) + Pm, O' - 1)] + If 2 ^ 2 }, 

p^(') = °) A ^' {p^(‘ " 0 + T [Pm, 0 _ J ) + Poxo(' ~ 1)1 + I t2<7 ^’ 


(26a) 

(26b) 

(26c) 

(26d) 

(26c) 

(26f) 

(27a) 

(27b) 

(27c) 


and 

It is clear that significant simplifications are required in order to find a steady-state solution. We will consider the 
limit that o w becomes much larger than all the other measurement and process noise variances. This means that the 



filter is ignorant of the spacecraft dynamics, in the sense of modeling the torques on the spacecraft and its response to 
these torques. In fact, we may have some knowledge of these quantities, but the filter does not make use of this 
knowledge, instead relying entirely on the RIG for information on the angular motion. This is the most common 
case in practice. Equation (27d) shows that Pim {i) becomes large along with a w . so we will take 


p wa> d) = aj^i + p^d) 


(28) 


where the tilde denotes the part that remains finite as cr„ becomes large. None of the other components of P(i r+) 
contain terms that become large in proportion to cr w . Now Eq. (25) gives 


A, =T 3 (a M +i)<T»+A i , 


(29) 


where A, is the finite part 

A, s p"(i-l) + 2r[p„(i-l) + p^(i-l)] + r 2 [p bb (i- l ) + 2 Pba,('- l ) + Po»,( i -W + Ta > + * r3<T « + <T ‘ ' 
Substituting Eq. (28) and 

A' 1 = r- J (a w + i)-'tr; 2 [ 1 - f" 3 (a,-i + i)" 1 ^ +...]. 
into Eqs. (26) and (27) and equating terms of first order in ff 2 gives 
a, =(i + a ,-,)/( 1 + 3 a i-i)- 

Equating terms that are independent of cr 2 gives the ten equations 

PeeiO = Peed ~ D + P» d - 1 )~ *P»V ~ « ‘ + 2 l W' -1 )“ />•<*' " ^ + !) + + ^ + ' 

P»{ i) = Petd - 1) - p»d - 1) - T Pu d - 1) - i^ 2 a 2 , 

P kb d) = P^d-D + ra], 

P<»d) = cr, 2 . 

Pb*d) - 0, 

P*d) - g ] • 

Pea d) = Ptod - !)■ - p«,d - !)• - TPtod - 1 )- + > *<-, r-'{p H d - D- P»d ~ D 

+ r[ Pfa d - 1) + 2 Pif d - 1) - Pem d - 1) - Petd - 1)] + f 2 [P^d-D + P^d - 1)1+ r<x 2 + 

Pbu d) = Pbvd - 1) - KAP»d - 1) + p b „d - 1) + r-' Pbf d - D + ^tr 2 ], 

and 

- Paxa (i) = (*_, - D 2 Pm d - 1) + 2* j .,(* I ,, - Off 1 P«,d - 1) + Ptcd - 1)1 

+ * 2 ,t‘ 2 [/j w (< - 1) + 2 rp M (( - 1) + - 1) + T(T 2 + + CT 2 1. 

where 

*, =(i + a,)/(i + a,>- 

Terms that go to zero as o w becomes infinitely large are ignored. 


(30) 

(31) 

(32) 

(33a) 

(33b) 

(33c) 

(33d) 

(33e) 

(33f) 

(34a) 

(34b) 

(34c) 

(34d) 

( 35 ) 


We will see below that the star tracker update doesn't affect a,. Thus we can assume that the coefficient a has 
attained a steady-state value resulting from repeated iterations of Eq. (32). This is the positive root of 

a = (| + a)/(l + 3a), ( 36 ) 


or 


a = 1 /V 12 . 

The corresponding steady-state value of k is, from Eq. (35), 

k = 3-V3. 

We now want to consider the overall effect on the covariance of the n gyro updates. Since P to , P*^ and p, 
do not appear on the right side of Eqs. (33), these equations can be combined in the 3x3 matrix form 


(37) 


(38) 
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(39) 


This equation represents the combined effect of propagation over a time step and a gyro update. It is easy to show by 
induction that for i > 0 

-|T 
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(40) 


tr 2 0 

Star Tracker Measurement Processing 

Using the notation and steady-state relation introduced at the beginning of the previous section, Eq. (40) gives 

P„(0-) = P„(flT+) ^ 2 (41a) 

= p„(0+) + P„(0+) - 2P*(0+) + 2 T [P w (0+) - P»(0+)] + T 2 P bb (0+) + To]+\ T\ T 2 + o ) . 

Pee (0-) = P» («T+) = Pee (0+) - P* (0+) ~ T P bb (0+)~i T 2 o ; , ( 41b ) 

P»(0-) = P bt ,(nr+) = />»( 0+) + Tct 2 , 

P^(0-)= ^OlT+) = CT 2 , 

P w (0-) = P*0iT+) = 0. 


(41c) 
(4 Id) 
(4 le) 


and 


P„(0~) = P„(nT+) = cr 1 ' ■ 


(4 If) 



( 42 ) 


The star tracker measurement update at time zero gives 

P( 0+) = P( 0-) - P(0-)Hl[H„P(0-)Hl + a 2 ]' 1 H s , P( 0-). 


Inserting this into Eqs. (41a— c) gives the steady-state conditions 

[a] +TP eb (0-)? = [P ee (0-) + o;][T 2 P bb (0-)-P ee (0-) + T(T; + \T'ol + 2 a;+o;] 

[a; - ^(0-)+ r/^O-)]^-) = [F*( 0 -) + a 2 ][7/> w (0-) + jr 2 <x 2 ]. 


(43a) 

(43b) 


and 

/ , i( 0 -) = [/ , w (0-) + cr„ 2 ]ra u 2 . (4 

The steady-state values of P ^ , />*, , and P^ are not of interest. We note, however, that Eq. (42) does not change 
the component of P^ proportional to cr 2 , so a is left unaltered by the star tracker measurement, as promised. 

Defining a dimensionless variable £ by the relation 


Eq. (43c) gives 

P»(0-) = (C 2 -l)cr 2 , 

and it follows from Eq. (42) that 

P„(. 0+) = (l-CVn- 

Equations (43b) and (42) give 

P»(0T) = [? ■ - (1 + S])C )T- ll2 o u a n ±\Ta] , 

where 


S t =aja„. 

We similarly define, following Farrenkopf, 

S^T^oJo, 


(45a) 


(45b) 

(46) 


(47a) 


(47b) 


and 


= T 3 / 2 / 


S =T 


J°*- 


The remaining task is to find £by solving Eq. (43a), which can be written 

[C 2 -2(|S,+y)C + l + 5, 2 ][C 2 -2(iS u -y)C + l + ^] = 0, 

with 

y = (l + S 2 +iS 2 +iS 2 ) ,/2 . 

The physically significant root of Eq. (48) is the largest root, 

C = y + iS„+!(2yS.+S 2 +*S 2 ) ,/2 . 

where yis assumed to be positive. Equation (50) gives 

(1 + s;)C' = r + is. - ±( 2 y + s 2 +4s y 2 ) l/2 . 


(47c) 


(48) 

(49) 


(50) 


(51) 



which can be used to rewrite Eq. (49) in the more convenient form 

P bb (0+)^(2yT' ,: o u o n + ct ; + a u ±iTcr u 


( 52 ) 


The principal results of this paper are given by Eqs. (45). (49). (50) and (52). Although the delation is much more 
involved than Farrenkopf s and the notation is somewhat different, the final results only differ by the 5, term in q. 
(49). For <7, = 0. these equations are completely equivalent to FarrenkopFs. It is remarkable that these results do not 

depend in any way on the gyro update interval r. 

Approximate Forms 

Approximate expressions valid for rapid updates are obtained by retaining only the lowest order terms in S v and S u . 
Equations (49) and (50) give 

(53) 

(54a) 
(54b) 


T/2 


? = (1 + s; y 1/2 + ± [2 (1 + s/ ) 1/2 S u + s v 2 ] 

which gives the approximate forms for the angle and drift bias variances 

P ee (0-) = (a 2 n +a 2 e ) ll2 {2T if2 a u (al+a 1 t ) l,2 + Tal] l/2 + a 2 t , 

Pee( 0 +)» 


a - 4 [2 7 ♦ T *]* + 


>.*+< t 2 ,) 312 


and 


P hb (0+) = [2T ]/2 o u (al+G z f ) 


2\l/2 


+ CTy] 


. 2 , 1 / 2 , 


(55) 


The difference between the pre-update and post-update drift bias variance is of order a u , which is negligible in this 
approximation. When (7, is negligible, the post-update and pre-update angle variances are equal in the approximation 

of rapid updates, also. 


Numerical Example 


Consider a rate- integrating gyro with 

cr v = 0.025 deg/ Vhour = 7.27 // rad/ Vsec, 

<7 = 3.7 x 10' 3 deg/hour 3/2 = 3 x 10" 4 /irad/sec 3/ \ 


(56a) 

(56b) 


and 


a f = 15/irad. 


(56c) 


These numbers are characteristic of a ring-laser gyro with very low drift but with significant angle white noise. 
Assume that the star tracker measurement noise is 

a n = 15 jUrad . 


(56d) 


Figure 1 shows the steady-state pre-update and post-update angle standard deviations, which are the square roots of 
/WOT), for star tracker update times between 0.01 sec and 100 sec. The solid curves are the exact values given by 
Eqs. (45), the dashed curves are the result of ignoring a, in these exact equations, and the dot-dash curves are rapid- 
update approximations of Eq. (54). For each pair of curves, the upper curve is the pre-update value, and thelower 
curve is the post-update value. The figure shows that the rapid-update approximations always lie between the p 
update and post-update exact values, and that these approximations are reasonably accurate for update times ot ess 
than one second with this set of gyro parameters. 


The steady-state drift bias standard deviations are not plotted, since they always lie between 0.0467 /tradians/sec and 
0.0468 /tradians/sec for these gyro parameters, whether or not the a, terms are included. The lower value is just 
equal to ^ct u ct v , the common limit of Eqs. (52) and (55) for small <7„ . 


State Update Equations 


The steady-state Kalman gains and state update equations are also of interest. The steady-state gains for the gyro 
update are given by 


K s = P(iT+)Hja' r 1 = [1, 0, 1, k/rf . 


(57) 


It is remarkable that these gains are constant, independent of the number of gyro updates since the previous star 
tracker update. The state update equation is 

x(ir+) = x(iT-) + K g l<t>i 0Of-)]. 



Figure 1. Steady-state pre-update and post-update angle standard deviations 

Solid curves from Eqs. (45). dashed curves from Eqs. (45) with cr, = 0. dot-dash curves from Eq. (54) 




Thu gives, for the individual components: 

00 r+) = d(iT-) + 0 I -0(ir-) = 6((i - l)r+) + 0, - <9(0 - 1K+) _ T ^((' “ ^ 59a ^ 

(59b) 

6(iT+) = 6(iT-) = 6((i-1)t+). 

1,- \ A (59c) 

0(ir+) = <90 r-) + <9, - 0 (<t-) = (9, . 

and 

^ (i r+) = dHiT-) + Arr" 1 [0, - ^(ir-)] = (1 - *) ai ((i - l)r +) + [0, - ^C(i - l)r+) - r * «i - l)r+)] , (59d) 

where we have used Eqs. (10) and (1 1) to include the state propagation from the previous gyro update. 

The first three of Eqs. (59) have a straightforward interpretation for i > 2. In this case <9(0 - 1 ) T+ ) = 0-i ^ the 
estimate of the spacecraft rotation angle is the previous estimate plus the rotation sensed by the gyro, corrected for 
the estimated drift bias. The estimate of the bias is not refined by the gyro measurement, as expected; and the 
estimate of the gyro-sensed angle is equal to the gyro measurement itself. We will return to the i - 1 case after 

considering the star tracker update. 

Equation (59d) is more difficult to understand. For k equal to unity, it would just say that the best estimate for the 
angular velocity is the gyro-sensed rate over the last measurement interval, corrected for drift. For k less than one, 
the updated angular velocity estimate would be some average of the previous estimate and the gyro-sensed rate over 
the last measurement interval. The steady-state k is greater than one, however, which means that the previous 
estimate of the angular velocity is given a negative weight in this average. It is worth noting that the variance of the 
angular velocity estimate is very large in our model, as shown by Eq. (28). 

The steady-state gains for the star tracker update are given by 

i-r 2 

,t_-2 -(CTT'S. (60) 


K„ = P(0+)//X 2 = 


(SJ O 2 ’ 


where we have used 


P»( 0 +) = r 2 P»( 0-) = -C-'T' l2 o u o a 


The state update equation is 


p*(o+> = r 2 v°-) = r 2 ®?- 


jc(0+) = i(0-) + K„ [0 O - 0(0-)]. 


Star tracker updates of d ) , which depend on the steady-state value of ^(0+), are not generally performed in 
practice. The other three components of the state vector are updated by 

0(0+) = 0(0-) + (1 - C 2 )[0 O - ®(0-)l. 

6(0+) = 6(0-) - (C7T 1 S„[0 O - 0(0-)] , 


0(0+) = 0(0—) + (5, /O 2 [®0 - ^(° _ )1 • 


(63c) 


It is convenient to define a modified angle estimate 0 by 

0(if+) = 0(/T+) 

for i = 1.2 «, and 

0(0+) 2 0(0+) + 0 O -0(0+). 


(64a) 


(64b) 


where 0„ IS understood to denote the gyro measurement immediately preceding the star tracker measurement at time 
zero. Thus 8 is the optimal angle estimate following a gyro update, but differs from the optimal estimate after a star 
tracker update. Now Eq. (59a) can be written 

0(zr+) = 0((« - 1 )t+) + 0, - </>,_, - r b((i - 1)t +) 

for all gyro updates. With this restructuring of the update equations, the need to estimate 0 vanishes. 

Since 0(0-) is the angle estimate following the gyro update immediately preceding the star tracker measurement, it 

is consistent with our labeling conventions to denote this quantity by 0(0-) on the nght side of Eqs. (63). This and 
Eqs. (50), (51), (63), and (64) enable us to write the star tracker update of the optimal angle estimate as a two-s ep 

process 

(66a) 


0(0+) = 0(0-) + r'(2y + s; + is u 2 ) ,/2 [0 o - 0(0-)], 


followed by 


(66b) 


0(0+) = (cr„ 2 + 0(0+) + o) 0 O ]. 

Equations (63) and (65) show that the update of Eq. (66a), which vanishes if <T„ = <T V = 0 , is propagated into future 
estimates of the angle and the drift bias. In contrast to this, the update of Eq. (66b), which does not depend on a or 
a and vanishes if a e = 0 , is only effective at the time of the star tracker measurement. The appearance of Eq. ( ) 

as the optimal combination of independent quantities 0(0+) and 0 O with standard deviations <x, and a. is very 
misleading, since these quantities are correlated by Eq. (66a). 

Conclusions 

Analytic expressions for the steady-state accuracy of a single-axis Kalman filter combining data from a gFo andan 
angle sensor have proven very useful in tailoring attitude sensor specifications to mission requirements. This paper 
very simple modification Famenkopfs analyuc expressions to include the effect of *£ “* 

rm the output of a rate-integrating gyro. The resulting equations depend on the 

independent of the gyro update interval. The effect of the angle white noise is significant for a specific set of gyro 
^aSSm Inalyzed. A ra^d-upda^ approximation to die stead, -state accuracy is reasonably accumte for sou hacker 
update rates of once per second or faster for these gyro parameters. 

The model developed in this paper leads to an effective two-component filter for the attitude angle and gyro drift bias, 
much like Farrenkopfs. In fact, the filter is identical with Farrenkopfs if gyro angle output noise is absent^ In the 
presence of gyro angle output white noise, the update of the angle estimate resulting from a star tracker measurement 
can be broken into two parts. The first part, which vanishes when the gyro drift parameters are zero, is propaga 
forward into future angle estimates. The second part of the update, which is independent of the gyro t parame ers 
“ishesTf the gyro readout white no.se is zero, is on.y effective at the time of the star tracker measurement. 
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